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We present a first-principles study of effects of small biaxial strain (|e| < 5%) and perpendicular 
electric field (E-field) on the electronic and phonon properties of low-buckled silicene and germanene. 

With an increase of the biaxial strain, the conduction bands at the high symmetric F and M points 
of the first Brillouin zone shift significantly towards the Fermi level in both silicene and germanene. 

In contrast, the E-field changes the band dispersions near the F and open a small band gap at the 
K point in silicene. We found that the field-induced gap opening in silicene could be enhanced by 
a compressive strain while mitigated by a tensile strain. This result highlights the tunability of the 
electronic structures of silicene by combining the mechanical strain and the electric field. 


I. INTRODUCTION 

The success of graphene has stimulated great interest 
in novel two-dimensional (2D) atomic crystals for inter¬ 
esting physics and diverse applications^. Silicene and ger¬ 
manene, the silicon and germanium analogs of graphene, 
are currently attracting growing attention^ - — . Similar to 
graphene, both of these crystals are shown to exhibit lin¬ 
ear band dispersions near the Fermi level at the K point 
of the first Brillouin zone (BZ) in the absence of spin- 
orbit coupling (SOC)£. Due to their low-buckled struc¬ 
ture, the sublattice inversion symmetry can be broken by 
applying a perpendicular electric field (E-field), leading 
to a sizable band gap of up to tens of meV— . Such 
a field-tunable band gap is useful for possible electronic 
device applications. The recent experimental fabrication 
of silicene on different substrates^ - — 7 has facilitated to 
further probe the material and exploit its properties. 

Since silicene and germanene are most likely to be fab¬ 
ricated on substrates^ - —, external strain induced by 
lattice constant mismatch between silicene/germanene 
and substrate may be present. In a recent theoretical 
proposal^, silicene may be grown on graphene, with a 
compressive strain up to 5 = —3% on the silicene lattice. 
Therefore, understanding of the strain effect on the elec¬ 
tronic and phonon properties will be essential for silicene- 
based application s 19 ! 20 . On the other hand, the in-plane 
strain provides an independent degree of freedom to tune 
the electronic properties in these 2D crystals^ 23 . In 
fact, strain has already become a powerful means of 
tailoring the electronic structure and to affect the car¬ 
rier mobility in silicon-based material s 24 ! 25 and in two- 
dimensional graphene^ - —. It has been reported that the 
in-plane biaxial and uniaxial strains could dramatically 
change the Fermi velocities in silicene^. Small biaxial 
strain of 5% also opens up channels possibly useful for 
an enhanced electron-phonon coupling in silicene^. Fur¬ 
thermore, strain may introduce delicate change on the 
phonon dispersions, dictating interesting Raman spectra 
as observed in graphene^! - — . 


It is expected that by combining an in-plane strain 
with a perpendicular E-field, a wider range of tunabil¬ 
ity on the electronic structure of silicene and germanene 
may be possible. How will these two factors affect each 
other? To address this question, we have carried out 
detailed first-principles calculations to study electronic 
properties of silicene and/or germanene with small biax¬ 
ial strain (with (|er| <5%)) under an additional E-field up 
to 5 V/nm. A possible E-field strength of up to 3 V/nm 
in bilayer graphene has been realized using dual-gated 
structure in experiment^. We limit our study to this 
small strain range so that silicene and germanene will not 
alter the semimetallic nature. Besides, such a small strain 
will not change the crystal structure dramatically, and 
achievable in experiment s 37 ! 38 . Despite of recent work 
on silicene under strain^ 23,2 *, a first-principles study 
including both strain and E-field in silicene has not been 
reported. 

The paper is organized as follows. Calculational de¬ 
tails are given in Section II. In Section III.A, we will 
first consider the strain effect on the crystal structure of 
silicene and germanene. The strain effects on the band 
structures and phonon dispersions of silicene and ger¬ 
manene are presented in III.B and III.C. We then focus 
our discussions on silicene with both strain and E-field 
in Section III.D. A summary is given in Section IV. 


II. METHODS 

Our first-principles calculations were performed us¬ 
ing density functional theory (DFT) and density- 
functional perturbation theory (DFPT)i£ with the lo¬ 
cal density approximation (LDA) as implemented in 
the Quantum ESPRESSO code^i. Norm-conserving 
pseudopotentials^ for Si and Ge were adopted to de¬ 
scribe the core-valence interactions. The wave functions 
of the valence electrons were expanded in plane waves 
with a kinetic-energy cutoff of 36 and 40 Ry for silicene 
and germanene, respectively. A Monkhorst-Pack uniform 
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FIG. 1. (Color online) (a) Top and (b) side views of the 
low-buckled silicene crystal lattice. The lattice constants of 
unstrained and strained silicene are denoted as ao and a , re¬ 
spectively. The buckling constant Ao is shown in (b). (c) The 
high symmetric points T, K, and M in the first BZ of silicene. 

/c-grid of 36 x 36 x 1 is used to do the self-consistent calcu¬ 
lations. For the phonon dispersions, we used a 6x6x1 q- 
grid. A vacuum region of 20 A is introduced to eliminate 
any artificial interaction between neighboring supercells 
along the perpendicular direction. The relaxed lattice 
constants are ao = 3.83 A for silicene and 3.95 A for 
germanene. The relaxed buckled separations are Ao = 
0.44 A for pristine silicene and 0.64 A for pristine ger¬ 
manene. These results are in good agreement with pre¬ 
vious DFT calculations^^. After we obtained the struc¬ 
ture of unstrained silicene and germanene, various biaxial 
strains were applied to the siliecene and germanene lat¬ 
tices, respectively, with the biaxial strain e defined as e 
= (a — ao)/ao- Here a and ao are the strained and un¬ 
strained lattice constants, respectively, as schematically 
indicated in Fig. [lj According to this definition, a nega¬ 
tive 5 means the compressive strain, while a positive value 
indicates the tensile stress. For each strain, the atomic 
positions were fully relaxed until the force is smaller than 
0.02 eV/A. 


III. RESULTS 

A. Effects of strain on the buckling constants 

In silicene and germanene, the low-buckled crystal 
structure is the main feature distinct from the planar 
graphene. This buckling is characterized by the vertical 
distance between the two silicon atoms in the unit cell, 
i.e., the buckling constant Ao, as indicated in Fig. [ljb). 
It is expected that an in-plane tensile strain decreases Ao. 
Figure [2] shows Ao/2 as a function of £ in both silicene 
and germanene for e =-0.01~ 0.05. When a positive bi¬ 
axial strain is applied, the lattice constant increases. As £ 
increases from 0 to 0.05, Ao/2 decreases 41% from 0.22 to 
0.12 A in silicene. In contrast, the buckling constant de¬ 
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FIG. 2. (Color online) Half of the buckling constant Ao/2 as 
a function of applied strain £ in silicene and germanene. 


creases only 6% in germanene for the same strain range, 
much smaller than in silicene. This is mainly due to the 
distinct atomic radii of Si and Ge. The single-bond radii 
are 1.09 A for silicon and 1.22 A for Ge, respectively. 
The larger atomic radius of Ge makes the Ge-Ge bond 
insensitive to a relatively small strain. These results are 
overall consistent with previous DFT calculations^. 


B. Effects of strain on the band structures 

The biaxial strain has dramatic effects on the electronic 
properties of 2D crystals. Figure [3] shows the band struc¬ 
tures of silicene and germanene as a function of strain 5 
from -0.01 to 0.05. In this range, a few interesting fea¬ 
tures can be summarized for silicene. First, the linear 
band dispersions near the K point keep almost intact for 
both silicene and germanene. Second, the biaxial strain 
shows more evident effects on the bands at the Y and 
M points than at the K point. In particular, the con¬ 
duction band minimum (CBM) moves towards the Fermi 
level as 5 increases. The CBM decreases from 2.4 eV for 
5 = -0.01 (Fig.[3fa)) to 0.5 eV for 5 = 0.05 (Fig.[3fd)) in 
silicene. This band also exhibits anisotropic dispersions 
along the Y-M and Y-K directions. In contrast, the va¬ 
lence band maximum (VBM) at Y moves downward into 
deeper energy by the tensile strain. In Fig.[4]we show the 
isosurface of the charge density for the CBM at Y in both 
silicene and germanene. Interestingly, this band exhibits 
somewhat 7r-like bonding features. 

The biaxial strain has significant effects on the 
parabolic conduction bands near the M point. The 
parabolic bands moves from 1.2 eV for s = -0.01 
(Fig. (3ja)) to 0.7 eV from 5 = 0.05 (Fig. (3^d)) in sil¬ 
icene. The parabolic band begins to cross with the i r* 
band near the M point, as shown in Fig. [3^d). 

Comparing to silicene, the strain shows similar but 
more evident effects on the band dispersions in ger¬ 
manene. The CBM at Y is very sensitive to the strain. 
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For 5 = 0.03, the CBM at T already touches the Fermi 
level, leading to semimetal-to-metal transition, as shown 
in Figs. [3] (g) and (h). Similar result has also been re¬ 
ported for silicene with larger tensile strain^. 

The strain-induced shift on the parabolic conduction 
bands (band No. 6) at the M point changes the band 
separations and consequently opens new scattering chan¬ 
nels for the electron-phonon coupling (EPC). It has been 
reported that there will be a significant enhancement 
of the EPC in electron-doped strained silicene^. To 
see whether the sing-particle band shift induced by the 
strain is still valid after including the many-body ef¬ 
fects, we have performed GW calculations using the 
ABINIT code^. Following the standard approach of one- 
shot GoWo method, Kohn-Sham eigenvalues and eigen¬ 
functions were firstly obtained by DFT-LDA calcula¬ 
tion and then used as a starting point to do the GW 
correction^ - —. Due to the non-local nature of the self¬ 
energy operator, the K-point used for the GW calculation 
has to belong to the K-point grid that has been chosen for 
the Kohn-Sham self-consistent calculation^. A 18 x 18 x 1 
Monkhorst-Pack grid is used. The GW band structures 
were interpolated based on the GW energies explicitly 
calculated at the grid and the DFT-LDA band disper¬ 
sions. The screening in the GW calculation is treated 
with the plasmon-pole mode l 44 ! 45 . The polarization func¬ 
tion is calculated within the random phase approxima¬ 
tion. Coulomb interaction is truncated using a cut-off 
radius of 5 A along the perpendicular direction. The 
number of bands used to calculate the screening and the 
self-energy in the GW method is chosen to be 138. The 
cut-off energy of the plane waves is set to 327 eV (24 Ry) 
to represent the independent particle susceptibility and 
381 eV (28 Ry) to represent the dielectric matrix and to 
generate the exchange part of the self-energy operator. 
The results have been shown in Fig. [3j 

From Fig. [3j one can see that overall the GW band 
structures follow the DFT-LDA trend, i.e., the parabolic 
bands at the M point shift towards the Fermi level as 
e increases. The GW quasiparticle bands also exhibit 
larger Fermi velocity near the K point, in agreement with 
previous report^. From Fig.[3j the band separations will 
become to match the phonon energies for even smaller 
tensile strain. Different from the single-particle picture^, 
one may expect that the EPC enhancement will occur at 
an even smaller tensile strain. 


C. Effects of strain on the phonon dispersions 

The effects of biaxial strain on the phonon dispersions 
in silicene and germanene are shown in Fig. [5] Note that 
in Figs. 5(a)-(b) and 5(e)-(f), small negative phonon fre¬ 
quencies (within -20 cm -1 ) appear on the acoustic branch 
along the Y-K direction. These small negative modes, 
however, do not mean the mechanical instability of the 
system. Instead, it mainly comes from the interpola¬ 
tion process as implemented in the Quantum Espresso 


code^i. Hence, both silicene and germanene will be sta¬ 
ble with the biaxial strain 6 from - 0.01 to 0.05. Generally, 
isotropic compression results in phonon stiffening of the 
vibrational mode, while isotropic tension results in a de¬ 
crease on the vibrational frequency (phonon softening). 
From Fig. [5j the tensile strain lowers the vibrational fre¬ 
quency, especially the in-plane and out-of-plane optical 
modes in silicene and germanene. 

The rate of change in frequency as a function of strain 
for a given phonon mode in a crystal is determined by 
its Griineisen parameter 7 ^. It is a crucial parameter to 
quantify the rate of the phonon mode softening (stiffen¬ 
ing) under tensile (compressive) strain and determines 
the thermo-mechanical properties. In metrology applica¬ 
tions, accurate determination of the Griineisen parame¬ 
ters are key for quantifying the amount of strain in the 
system. In presence of biaxial strain, the Griineisen pa¬ 
rameter 7 g for a particular band m associated with an 
in-plane phonon mode read s 31 ! 35 ! 47 : 


_ 1 dujffi 

19 ~ ~2< We 


( 1 ) 


where e is the biaxial strain applied to the system, and 
and uj m correspond to the phonon frequencies at zero 
strain and in presence of an applied strain, respectively. 

We have performed a linear fitting of the phonon fre¬ 
quency shifts as a function of 6. The Griineisen param¬ 
eters are then calculated using Eq. (1) based on the ob¬ 
tained linear slopes. Table I summarizes the results. For 
comparison, we also listed the calculational results for 
graphene. The calculated 7 g for the Y — E^ g and K -A[ 
modes in graphene are 1.8 and 2.8, in good agreement 
with the values of 1.8 and 2.7 obtained by Mohiuddin 
et a/4I. Surprisingly, silicene has the smallest 7^ of 1.3 
for the T — Eg mode, which is 14% smaller than in ger¬ 
manene and 25% smaller than in graphene. Our LDA 
value is smaller than the GGA result of 1.6 obtained by 
Kaloni et a/. 23 . As can be seen from Table I, the 7 ^ for 
the K -A mode are nearly 2.8, 2.0 and 2.1 for graphene, 
silicene, and germanene, respectively. We conclude that 
the 2D band in the Raman spectra, i.e., the overtone of 
the D mode (corresponding to the K -A mode), may be 
more sensitive to the biaxial strain than the G band. 

Finally, we consider the effects of strain on the EPC 
for the T — Eg mode in silicene. In our previous work^S, 
the EPC in germanene is found to be an order of mag¬ 
nitude smaller than in silicene. Therefore, we only focus 
on silicene and compare the results with graphene. Since 
at k = K, the electronic states are doubly degenerate, 
we calculate the average EPC square over the Fermi sur- 
face defined as {g 2 q ) F = JV • \g( K + q )i,Kj\ 2 / ^ with q = T, 
where the sum is performed over the doubly degenerate 
7 r bands at Ep. This quantity is found to be a good 
indication of the EPC strength of the E g mode for the 
electronic states near the BZ corner. In silicene, the cal¬ 
culated value is (g^F = Yli j \dKi,Kj\ 2 / 4 = 0.0223 eV 2 . 
As depicted in Fig. \6\ we found that the EPC square 
decreases as 5 increases. At e = 0.05, the EPC square 
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Silicene Germanene 



FIG. 3. Band structures of silicene (left column) and germanene (right column) as a function of strain. The Fermi level has 
been shifted to zero. Solid line: DFT-LDA band dispersions; Dotted line with symbols: interpolated GW band structures. 



TABLE I. Calculated Griineisen parameters for the highest 
optical modes at Y and at K in graphene, silicene and ger¬ 
manene, respectively. 


Graphene 

Silicene 

Germanene 

r - E 2g K-A 1 ! 

r — Eg K -A 

Y-E g K-A 

1.75, 1.8^ 2.75, 2.7 a 

1.31, 1.64^, 1.95 

1.52 2.11 


a Ref. HI]. 
b Ref. 0. 


FIG. 4. (Color online) Charge density distribution for the 
CBM at T with strain e = 0.05. (a) Isosurface and (b) cross- 
sectional view of the charge distribution in silicene; (c) Isosur¬ 
face and (d) cross-sectional view of the charge density distri¬ 
bution in germanene. The cross-section planes are indicated 
by dashed line in (a) and (c), respectively. 


in silicene is 0.0197 eV 2 , about 12% smaller than in un¬ 
strained silicene. In graphene, the EPC square is 0.0409 
eV 2 , nearly 7% smaller than the unstrained one. This re¬ 
sult shows that the EPC for the Y — E g mode in silicene 
is almost two times more sensitive to the biaxial strain 


than in graphene. 


D. Effects of both strain and E-field on the band 
structures 

The perpendicular E-field has been shown to open a 
small energy gap in silicene££. This feature makes sil¬ 
icene a very interesting platform for electrically tunable 
device applications. However, the required E-field is too 
strong (up to 5 V/nm), and the opened band gap is rela¬ 
tively small (~ 20 meV), which limits the practical appli¬ 
cation of silicene. Will the biaxial strain possibly enhance 
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FIG. 5. Phonon dispersions of silicene and germanene as a 
function of strain e. Left column: silicene; Right column: 
germanene. 
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FIG. 6. EPC square as a function of strain in silicene and 
graphene. 


-graphene 

-silicene 


-•-• 


-0.01 0.00 0.01 0.02 0.03 0.04 0.05 

Strain e 


the gap opening? Since the tunable strain range is rel¬ 
atively small in germanene, in this part we will mainly 
focus on the effects of both the E-field and the biaxial 
strain on the electronic structure of silicene. 

Figure [7| presented the band dispersions of silicene un¬ 
der various E-field (E = 1, 3, and 5 V/nm) for three typ¬ 
ical strains (s = -0.03, +0.03 and +0.05). When silicene 


is under a given strain 6, the additional perpendicular E- 
field shows two main effects: (1) The E-field moves the 
CBM at r down towards the Fermi level; and (2) the 
E-field opens a small band gap at K. The CBM at T will 
approach the Fermi level and eventually a semimetal-to- 
metal transition takes place at a high enough E-field (see 
Fig. [7]^a)-(c)). Meanwhile, the energy band gap opened 
by the E-field also increases as the field strength becomes 
stronger, as will be shown below. 

The effects of strain and E-field can be seen more 
clearly from Fig. [8j where the evolutions of the CBM 
and VBM at the T point are plotted as a function of 
the E-field for various strain s. Interestingly, the VBM 
doesn’t show any evident dependence on the E-field until 
5 V/nm. In contrast, the biaxial strain introduces rigid 
shifts on the VBM. In particular, a compressive strain 
shifts the VBM towards the Fermi level. 

On the other hand, the CBM at T exhibits a clear 
dependence on both strain and the E-field, as shown in 
Fig. [71 For a given E-field, the CBM shifts to higher 
energy from e = -0.05 to -0.01 and then moves back from 
5 = -0.01 to 0.05. For silicene with a specific strain £, 
the increase of the E-field will push the CBM downwards 
and touch the Fermi level at E = 5 V/nm. Moreover, 
the compressive strain seems to further enhance such a 
shift at T. 

In addition to the T point, the co-application of both 
the E-field and the strain in silicene has dramatic effects 
on the band gap opening at K. Figure [9] summarizes the 
results. 

The energy gap E g at K changes as a function of both 
the E-field {E = 0 ~ 5 V/nm) and the biaxial strain e 
(from -0.05 to 0.05), as shown in Fig. [9^a). For a given 
strain, the E g at K increases almost linearly as the E-field 
increases. For example, for 5 = 0.01, the E g increases 
from 0 to 30 meV when E changes from 0 to 5 V/nm. 
In contrast, the increase of the E g becomes saturated 
quickly as e changes to tensile strain. On the other hand, 
the compressive strain significantly enhances the band 
gap opening. For 5 = -0.05, the band gap changes from 
0 to 42 meV at E = 5 V/nm, more than 35% larger than 
that of 5 = 0.01—. 

The above behavior may be understood from the strain 
effect on the buckling constant Ao, as shown in Fig. [2j 
The compressive strain enlarges the vertical distance of 
two silicon atoms, and therefore, increases the on-site en¬ 
ergy difference between two sublattices. Using two-band 
model at the K point, one can write out the Hamiltonian: 


( +ei?Ao/2 0 \ 

V 0 -eUAo/2 ) ’ 


( 2 ) 


with E the field strength and Ao the buckling constant. 
The energy band gap is given by E g = eE Ao. From 
Fig. m the compressive (tensile) strain increases (de¬ 
creases) the Ao, and consequently, enhances (mitigates) 
the gap opening. As field strength increases, the screen¬ 
ing becomes stronger, especially for silicene under tensile 
strain with smaller buckling constant. This will decrease 
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(h) and 5 V/nm (i). 

FIG. 9. (Color online) (a) Energy band gap E g at the K point 
as a function of the E-field for various strain e in silicene. (b) 
The calculated enhancement factor a as a function of strain 
£. 


the displacement field strength and as a result, the band 
gap becomes saturated quickly, as shown in Fig. [9l This 
is similar to the results observed in bilayer graphene^. 

One can quantify the effect of strain on the gap open¬ 
ing. Since E g changes linearly with E, we define the 
enhancement factor a by fitting to the expression E g = 
a(e)E for each e. Here a(e) is a function of 5. E is 
the field strength in V/nm, and E g the energy band gap 
in meV. In Fig. [9fb), the calculated a is plotted as a 
function of 6. Clearly, a decreases almost linearly with 
respect to the biaxial strain e. This relationship can be 
described as a = 6.2-47.8^. For the compressive strain 
of 5 = -0.05, we have =1.4, while for the tensile 

strain of e = 0.05, the ratio is = 0.6. Clearly, for a 

given E-field, the compressive strain will enhance the gap 
opening by increasing a, while a tensile strain mitigates 
this trend. This result highlights the impact of the strain 
on the field-induced gap opening in silicene. 


IV. SUMMARY 


In summary, we have performed detailed first- 
principles calculations to study the effects of the biaxial 
strain and the E-field on the electronic and phonon prop¬ 
erties of low-buckled silicene and germanene. The small 
biaxial strain is found to dramatically change the conduc¬ 
tion bands at T and M, while the E-field mainly affects 
the bands at the T and K points. The Griineisen param¬ 
eters for the T — E g and the K -A mode in silicene and 
germanene are calculated and compared with graphene. 
We also show that the EPC matrix square for the T — E g 
mode in silicene is more sensitive to the biaxial strain 
than in graphene. Finally, the field-induced band gap 
is found to be sensitive to the strain. In particular, the 
compressive strain is able to significantly enhance the 
field-induced gap opening in silicene. 
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